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Abstract 

We present an ab-initio study of electron mobility and electron-phonon coupling in chemically 
modified graphene, considering fluorinated and hydrogenated graphene at different percentage 
coverage. Hexagonal Boron Carbon Nitrogen (h-BCN) is also investigated due the increased interest 
shown by the research community towards this material. In particular, the Deformation Potentials 
are computed by means of Density Functional Theory (DFT), while the carrier mobility is obtained 
according to the Takagi model (S. Takagi, A. Toriumi, and H. Tango, IEEE Trans. Electr. Dev. 
41, 2363 (1994)). We will show that graphene with a reduced degree of hydrogenation can compete, 
in terms of mobility, with silicon technology. 

PACS numbers: 73.63.-b,73.50.Dn,72.80.Vp,63.22.-m 



1 



Understanding electron mobility in chemically modified graphene is a matter of primary 
importance, in order to clarify if graphene will manage to come after industry requirements 
in the mid-term. 

Despite two dimensional graphene is now at the core of both fundamental and applied 



research, it presents lights and s 



radows, like its outstanding electronic properties, such as 



its extremely high mobility |2|, but, at the same time, the lack of an energy bandgap, 
which prevents its use in digital applications. 

Chemical modification of graphene has recently entered the research scenario as a possible 
candidate to tackle such a problem. Elias et al. 3| managed to experimentally induce 



an energy gap in hydrogenated graphene, previously predicted by Sofo et al. [4] through 
numerical calculations, while graphene fluorination 5j has shown a clear dependence of the 
band gap on fluorine concentration. 

An alternative option to obtain two-dimensional materials with tunable electronic prop- 
erties has been recently demonstrated by Ci et al. [(J, exploiting hexagonal Boron Ni- 
tride (h-BN) and including different percentages of C atoms (h-BCN), or (BN) n (C2) m 
(n, m = 1, 2 . . .) isomers. 

|— I 

However, fabrication techniques are still at an embryonic stage [3] , so that an experimental 
investigation of the main figure of merit for electronic application like electron mobility is 
still far from being performed. From this point of view computational approach is the only 
tool able to provide relevant information regarding the intrinsic mobility to be expected. 

To this purpose, in this work we focus on the electron-phonon interaction, in order to in- 
vestigate the ultimate intrinsic mobility /i in hydrogenated/fluorinated graphene at different 
percentage coverage, as well as in h-BCN structures. 

We follow the Bardeen and Shockley's [7] approach, where the atomic displacement asso- 
ciated with a long- wavelength acoustic phonon leads to a deformation of the crystal 8] , and 
in turn, to a shift of the electronic energy dispersion. The band edge shift is related to the 
differential displacement of the crystal by the electron-phonon coupling Hamiltonian js il]: 



H. 



el—ph 



£iV ■ u(r,t) 



(1) 



where E\ is the so-called deformation potential and u (r , t) is the displacement at the 
spatial coordinate f [12|. In order to evaluate V • -u(r, t), we define a second-rank tensor 
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dij = = e *i + fiji w here Ui and Xi are the i — th component of vector u and the cartesian 
coordinates, respectively, while and /y are the antisymmetric and symmetric tensors, 
respectively. In particular, describes a rigid rotation of the crystal and does not lead to 
a shift of the electronic energies, so that, to our purposes, it can be ignored, instead 
describes a strain induced in the crystal by the atomic displacements, which indeed is able 
to shift the bands. By definition, its trace is equal to the fractional volume change ^ of the 
elementary cell induced by the strain. Representing the long wavelength acoustic vibration 
by a plane wave with frequency u and wave vector q as u(r,t) = Suq sin (q ■ r — uit), the 
strain tensor fij in proximity of the T point, for the longitudinal acoustic (LA) mode, is 
diagonal, i.e. fa = qi5uoi, where qi and du^i are the i — th component of the wave vector 
q and the atomic displacement from equilibrium geometry, respectively. In this hypothesis, 
since q - u = we obtain V ;_u ~ ,n 
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For non-degenerate bands [13|, Eq.([T]) can be eventually reduced to 



5E (A:) = E x (2) 

where 5E (k) is the induced band edge shift due to acoustic phonon. Once known E\, 
the carrier mobility can be estimated by means of the Takagi formula [ijj], which reads 

where p is the mass density, T is the temperature, ks the Boltzmann's constant, e is the 

elementary charge, Si is the sound velocity, while m e is the mass along the transport direction 

(either m x or m y along the x and y direction, respectively) and is the equivalent density- 

of-state mass defined as rrid = y/m x m y . Si, m e and m<i are computed by means of ab-initio 

calculation, and are extracted from the phonon u(k) and electronic dispersion relations E(k), 

-l 



according to Si = 9uJ S^ and m e = h ° ^, ( 2 K> . The deformation potential E\ is instead 



dk 



d 2 E(k) 



dk 2 



computed through mimicking the lattice deformation due to phonons by multiplying the 
lattice constant by different factors (0.99, 0.995, 1.005 and 1.01), i.e. dilating and relaxing 
the primitive cell. Ei is then obtained computing the energy shift of the bottom of the 
conduction band 5E and applying Eq.(T5]). 



Ab-initio calculations have been performed by means of the Quantum Espresso 15] code 



using a plane wave basis set in the generalized gradient approximation (GGA) with the 
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Perdew-Burke-Ernzerhof (PBE) exchange correlation functional. A 50 Ry wave function 
cutoff and 400 Ry charge density cutoff have been considered, while the Brillouin zone has 
been sampled using a 30 x 30 x 2 Monkhorst— Pack grid. A 30 bohr layer of vacuum is 
considered to separate the sheet from its periodical images, which we have verified to be 
sufficient to avoid any unphysical interactions. The dynamical properties are calculated 
within the density-functional perturbation theory [161 ] . We have also verified that extracting 
the main parameter of interest, i.e. the effective mass, from DFT-GGA simulations manages 
to provide results in close agreements with those extracted from GW calculations (within 
5% of accuracy), but with fewer computational requirements (an overall simulation walltime 
almost 3 order of magnitude smaller). 

In Fig. [TJ we show the considered structures. Relaxed structures for 100% hydrogenated 
and fluorinated graphene and 50% fluorinated graphene (from here on H100%, F100% and 
F50%, respectively) are in agreement with results already shown in the literature 0, \v\ . 
25% and 75% hydrogenated graphene (H25% and H75%) have been studied considering a 
cell composed by 8 C atoms and 2 H and 6 H, respectively. In particular, the topological 
shown in FigJTJis the most energetically favourable over all possible geometric and spin con- 
figurations. Simple chemical consideration can support these results. The relaxed H25% 
configuration is the only one possessing a resonant chain of alternating double and single 
bonds, which guarantees structure stability. On the other hand, the considered H75% struc- 
ture shows a double bond inside the cell, which assures a stable structure. All the structures 
are nonmagnetic, since unsaturated C atoms are all neighbors and their p z orbitals form 
7r-bonding, which quenches magnetism. 

All the considered materials (except F50%, which is a metal) are semiconducting and 
have non-degenerate conduction band minimum at the T point. 

The effective masses shown in Tab. [I] are evaluated through a five-point second derivative, 
considering a k-point spacing smaller than O.OlA -1 in order to avoid nonparabolic effects. 
The sound velocities are extracted from the slope of the acoustic longitudinal phonon branch. 
H25% and H75% show different values for the effective mass along the YK and YM direc- 
tions, leading to two different values of electron mobility. The relative effective mass is quite 
large and close to one for the larger degree of functionalization. At the lowest degree of H 
coverage, the electron effective mass is small, since in the limit of no functionalization, the 
null mass of graphene has to be recovered. 
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FIG. 1: Geometric structures of the considered 2D materials: graphene sheet with different hy- 
drogen and fluorine coverage (a) 100%; (b) 75%; (c) 50%; (d) 25%; (e) BN and (f) BC 2 N. The 
primitive cell and the lattice vectors are also shown for each material. 

In order to understand the deformation potential E\ trend as a function of H percentage, 
we focus on the LUMO (Lowest Unoccupied Molecular Orbital) at the V point. 

In Fig. |2l we show the LUMO for the F100% and h-BN, corresponding to the largest and 
smallest obtained E\, respectively (the lighter the colour, the higher the probability density). 
In particular, the h-BN shows the LUMO strongly localized on single atoms (Fig. [2b), 
while in the F100% case the LUMO is spread over groups of atoms (Fig. EK). Since the 
deformation of the crystalline lattice due to longitudinal acoustic phonon is represented as 
an homogeneous, symmetrical dilation on the plane (picture at the top, where the direction 
of the stretching is depicted by the arrows in white), the LUMO of F100% is largely deformed 
by the phonon motion, leading to a stronger electron-phonon coupling, and consequently, to 
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TABLE I: Density mass p, sound velocity Si, relative effective mass m e , deformation potential E\ 
for the considered structures. 





P 


Si 


effective 


effective 


Ex 




( xW- 7 ^) 

\ J 


(xlO 4 ^) 


reduced 


reduced 


(eV) 








mass (TK) 


mass (TM) 




H 100% 


7.73 


1.78 


1.02 


1.02 


6.9 


H 75% 


7.65 


2.21 


1.59 


1.39 


8.58 


H 25% 


7.35 


2.46 


0.12 


0.53 


9.44 


v 1 nn% 

r luu/o 


1 8 8^ 


1.18 


0.53 


n ^ 

u.oo 




F 50% 


12.98 


1.34 


0.99 


0.99 


10.94 


h-BN 


7.59 


1.94 


0.97 


0.97 


3.66 






(TX) (TY) (TX) 


(TY) 




BC 2 N 


7.46 


2.1; 1.89 


0.39 


0.27 


18.74 



large value of E\. On the other hand, when the LUMO is strongly localized, the electron- 
phonon coupling is poor and the deformation potential E\ is consequently low, as in the 
case of h-BN. 

In Tab. [Til the electron mobilities computed for all the considered materials along the 
main symmetry directions are shown. H25%, H75% and BC2N exhibit different mobilities 
along the two considered directions, as a consequence of the asymmetry of the conduction 
energy band. 

The low mobility values for highly hydrogenated graphene are in quantitative agreement 
with available experimental data recently reported in the literature as a function of the 
exposure time to H |l8( . While fluorinated compounds show the smallest mobilities among 
the considered structures, H25% possess mobility comparable to bulk silicon, and an energy 
gap suitable for electronic applications, which could lead to a possible exploitation of such 
material in next-generation nanoscale devices. Finally, while BC2N shows poor mobility in 
both TX and FY directions, fi in BN is large enough for digital applications. 

To conclude, we have computed the electron mobility of hydrogenated and fluorinated 
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(a)100% fluorinated graphene 






FIG. 2: T-point LUMO wave function for (a) 100% fluorinated and (b) h-BN. In the picture at the 
top the direction of the stretching is depicted (white arrows). 



TABLE II: Electronic mobility [i of the considered materials along the main directions of symmetry. 

(TK) (TM) 

H 100% 104.8 104.7 

H 75% 45.7 52.1 

H 25% 3597.6 805.1 



F 100% 44.9 44.9 
F 50% 42.6 42.6 



BN 486.9 486.9 
(FX) (TY) 
BC 2 N 153.3 179.8 
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graphene as well as h-BCN from first-principles. The adopted method manages to provide 
accurate results with reduced computational requirements, since it avoids explicit calculation 
of electron-phonon scattering coefficients. Deformation potentials have been extracted from 
DFT calculations, which can be useful for future works dealing with Monte Carlo simulations 
of transport in functionalized graphene based transistors. Among the considered structures, 
H25% has better performance in terms of mobility, showing fj, comparable with bulk Silicon, 
which could open its exploitation as device channel in the next technological nodes. As a 
word of caution, we have to point out that the obtained mobility represents an upper limit for 
the intrinsic mobility achievable in such structures, since additional degree of inhomogeneity 
(not considered in the present work) could lead to mobility reduction. 
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216171), GRAND (n. 215752) grants, and by the MIUR-PRIN project GRANFET (Prot. 
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